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THE SOUTHERLY AND EASTERLY DEVIATIONS OF FALLING 

BODIES FOR AN UNSYMMETRICAL GRAVITATIONAL 

FIELD OF FORCE* 

BT 

WILLIAM H. ROEVER 

Introduction and Statement of Results. 

The adequacy with which theory may be made to account for a given class 
of physical phenomena depends largely upon the nature of the assumptions 
which are made in the mathematical formulation of the corresponding physical 
problem. This is well illustrated in the problem of the southerly deviation of 
falling bodies. In the treatment of this problem (and that of the easterly 
deviation) by Gauss, the assumption is made that the force in the statical 
field of force which is at rest with respect to the earth, i. e., the field of force in 
which the plumb-line is in equilibrium, is constant in magnitude and direction 
in the neighborhood of the path of the falling body. Under this assumption,! 
it follows from the laws of dynamics that the southerly deviation (denoted by 
S. D.) is given by the formula:} 

1 ¥ 

(1) S.D. = ^co2sin2<^— , 

<J 9 

in which w is the angular velocity of the earth's rotation, h is the height through 
which the body falls, and g and <^ are the acceleration due to weight and the 
astronomical latitude, respectively, at the place of observation. In a recent 
paper by the author § it is shown that, under the assumption of a distribution 
of revolution (i. e., that the potential function of the earth's gravitational 
(or weight) field of force is of the form f (r , z) , where r is the distance of a 
general point from the earth's axis of rotation, and z is that from a fixed plane 
perpendicular to the axis), the same theory yields for the southerly deviation 

* Presented to the Society December 2, 1911 (Southwestern Section) and April 5, 1912 
(Chicago Section). 

t Other assumptions are also made. The body is assumed to be a particle, and the effects 
of air currents, air resistance and the actions of the moon and sun are disregarded, 

t See Gauss's Werke, vol. 5 (1867), p. 502. 

§ These Transactions, vol. 12 (1911), pp. 335-353. 

Trans. Amer. Math. Soc. 30 469 



470 

the formula 
(2) 



w. H. roever: deviations 



S.D. = i[2co^sin2, + 5(0)J^^, 



[October 



in which dg / dx is the derivative of g along the meridian measured to the north, 
the subscript ( o ) denoting the particular value at the initial point Po of the 
falling body. It will be observed that when g is constant, dg j dx = , and 
therefore formula (1) is a special case of formula (2). 

For the potential function for which the Besselian ellipsoid is a level surface 
and the formula of Helmert gives the acceleration due to weight on this 
ellipsoid, the following table gives a few of the corresponding values of <i>, 
g , and dg / dx .* 

Table I. 






9 


^'< 


40° 
45° 
50° 
55° 
60° 


980.1457 
980.5966 
981.0475 
981.4847 
981.8949 


8.0399 
8.1568 
8.0259 
7.6517 
7.0463 



On the other hand, 

0)2 = 5.3173 • 10-» . 

The values given by formulas (1) and (2) are proportional to 
2 0)2 sin 2<i> and 2a)2 sm2<i> -\- 5{dg I dx)^. 
Since sin 2<^ never exceeds unity, it is evident from Table I, that, for the 




potential function for which the standard ellipsoid is a level surface, formula (2) 
gives values for the southerly deviation which are nearly five times as great 
as those given by formula (1). 

The expression {dgl dx)o/ go is the curvature, at Po, of the line of force 
which passes through Po of the statical field of force which is at rest with respect 
to the earth. Thus we see that the curvature of these lines of force influences 

* See 2nd reference in the second following footnote. 
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more strongly the southerly deviation of falling bodies than does the earth's 
rotation. 

Simple considerations show that the introduction of an attracting mass into 
a field of constant force, or into one whose lines of force have slight curvature, 
increases considerably the curvature of the lines of force in the neighborhood 
of the introduced mass.* 

This fact suggests that the proximity, to an experiment station, of a moun- 
tain or of a mineral deposit can change considerably the curvature of the 
lines of force of the statical field which is at rest with respect to the earth, 
and thus affect considerably the southerly deviation of falling bodies. But 
the proximity of such a disturbing factor would make the distribution of 
matter which produces the earth's gravitational field of force cease to be one of 
revolution. Therefore it becomes necessary to derive a formula for the south- 
erly deviation — one for the easterly deviation will also be derived — under the 
assumption of a distribution not of revolution. This has been done in this 
paper. (See formulas III, where |^ = S. D. and ^ = E. D.) It turns out 
that the first term of the expression for the southerly deviation, under this 
assumption, is the same as expression (2), and that the term in h^'^ , which 
does not appear when the assumption of a distribution of revolution is made, 
is negligible in comparison with the first term. 

The Hungarian physicist Baron Roland Eotvos has recently devised a 
method by means of which he can determine experimentally the derivative 
dgj dx, and also other second derivatives of the potential function W of the 
force due to weight.f The following table gives the observed values of these 
derivatives at a few stations which are about twenty kilometers east of Arad 
in Hungary and just west of some mountains. 



Table II. 



Station number 


10»^'^ 
dxdz 


<1 
dydz 


(aw aw\ 

'" V dy' ax^ J 


10' fT 

ax ay 


1,018 
1,032 
1,035 
2,159 
2,188 
2,192 


+32.9 
+25.8 
+15.5 
-37.0 

+77.8 
+82.3 


+88.1 
+62.8 
+69.8 
+34.5 
+59.8 
+64.5 


+37.0 
+72.8 
+25.5 
+34.3 
+21.5 
+28.2 


+ 3.8 
+ 2.5 
- 4.6 
-13.5 
+30.3 
+26.7 



* The opposite figure represents the lines of force of the gravitational field which is produced 
by the introduction of a mass m, at 0, into a field which originally was constant (and of which 
the lines of force were straight lines parallel to OY). The line AIB is the locus along which 
the curvature has not been changed; in this case, it is the locus of the points of inflection of the 
lines of force of the new field of force. 

tSee Encyklopadie der mathematischen Wissenschaften, Band VIib, 
§ 23, p. 166. Also, Verhandlungen der funfzehnten allgemeinen Con- 
ferenz der interna tionalen Erdmessung, 1 (1906), pp. 337-395. 
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Here x is measured to the north, y to the east, and z downward along the 
vertical. Hence 

dW dg d^W dg d'W 



9 = 



dz ' dx dxdz ' dy dydz ' 



For the latitude (46° 10') of this region, the normal values (i. e., the values 
corresponding to the potential function for which the standard ellipsoid is 
a level surface) of the quantities given in Table II are 

8.1, 0, 4.9, 0, 
respectively. 

At station 2,192 the value of the expression 2w^ sin 2<^ + 5 {dg/ dx)o is 
about 422 • lO"* , and hence formula (2) gives a value forty times as great as 
formula (1). On the other hand, at station 2,159, 2w^sin2(> + 5 {dg/ dx)o 
is — 174 • 10~*, and hence formula (2) gives a northerly deviation which is 
more than sixteen times as great as the southerly deviation given by formula (1). 
The local value of the derivative dg / dx = d^ W / dxdz thus appears to de- 
termine the extent and the sign of the southerly deviation of falling bodies. 

The easterly deviation of falling bodies is also influenced by the proximity 
of mountains or of mineral deposits. In fact, such disturbing factors yield 
the second term of the expression for the easterly deviation (denoted by E. D.). 
See formulas III. 



2,/^ h'" 5 {dg /dy)o 
3l/2a,cos,^^,-g ^^ 



(3) L. D. = 5 k2&) cos<^-[to— ^ « 



From Table II and the form of the second term in formula (3), it is evident 
that the effect of such disturbing factors is comparable in magnitude to the 
southerly deviation, but small in comparison with the first term of formula (3). 
The effect is also comparable in magnitude to the effect which Gauss finds is 
due to air resistance. Thus, for the data of Benzenberg's experiment in St. 
Michael's Tower, namely <i> = 53° 33' and h = 235 feet, the formula of Gauss 
for the easterly deviation (which is equivalent to the first term of formula (3)) 
gives the value 3.91 lines.* When air resistance is taken into account. Gauss 
finds for the same data the deviation 3.86 lines. Hence the effect of air re- 
sistance is —.05 lines. On the other hand the value of the second term of 
formula (3) for the data of station 1018 (Table II) is — .02 lines. Gauss 
states that for the southerly deviation the effect of air resistance is negligible. 
In this paper the effects of air resistance, air currents, and the actions of 
the moon and sun have not been taken into consideration.! While the com- 

* 144 lines = 1 foot. 

t Another matter which has not been taken into consideration in this paper is the weight 
of the string which supports the plumb-bob. In a field of force in which the lines of force 
are not rectilinear, this string will have some curvature. Therefore the position of a plumb-bob 
which is supported by a string with weight will be slightly different from that of a bob which 
is supported by a weightless string of the same length as the heavy string. 
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bined effects of these influences may be appreciable, it seems desirable to do 
first what has been done in this paper, namely to determine the effect of a 
local irregularity in the earth's gravitational field of force. The surprising 
conclusion is reached that known local irregularities in the earth's gravitational 
field of force (caused by the presence of mountains or of mineral deposits or even 
by large buildings or tunnels) (1) influence the southerly deviation to the extent 
of from — 16 times to + 40 times * the amount which is given by the formula 
of Gav^s for the southerly deviation {formula (1)), (2) affect the easterly deviation 
by amounts which are comparable with the effect which Gauss finds is due to air 
resistance. 

§ 1. Definitions of, and formulas for, the easterly and southerly deviations of 

falling bodies. 

In experiments for the determination of the easterly and southerly deviations 
of falling bodies, a plumb-line Pq R is supported at the point Pq from which 
a spherical body is later permitted to fall. On the horizontal plane which 
passes through the plumb-bob R, the easterly and southerly directions RE^ 
and RSg are drawn (Fig. 3). The plumb-line RPq is the vertical at ii.f 
The falling body, after being released from the position Po , moves in a path 
which, with respect to the axes R-E^, Sg, Pq, is the curve c . This curve 
pierces the horizontal plane E^RSji in the point C . The distance by which C 
lies to the south of RE^ is called the southerly deviation of the falling body, and 
that by which it lies to the east of RS^ is called the easterly deviation. 

In order to get expressions for these deviations the curve c is referred to the 
axes Po-v, I ) ? which pass through Po and are parallel to the axes 
R- Ejt, Sji, Po.X The equations of the curve c are (by Eqs. (17) § 7) 

5 



+ ^0, (sin ,>TF?f + cos <^ ( Wli + W^ - Aw')) f } , 



- W\ { ^ (4co2 sin <^ cos <^ + bW\i) t" 

-H ^ CO (9 sin <i,Wli + cos 4>Wl) f | , 



2 



* In the region of the Alps this range of values would probably be much greater. See note 
by M. Brillouin on the curvature of the geoid in the Simplon tunnel, Comptes rendu s, 
vol. 102 (1906), p. 916. 

t It is assumed that the plumb-bob is a heavy particle and that the line is weightless and 
perfectly flexible. 

t The axes Po — '^ , f , f do not coincide with the axes Po — t; , J , f . The latter axes 
are the easterly, southerly and vertical directions at Po . 



474 W. H. ROEVER : DEVIATIONS [October 

where t represents the time measured from the instant the body begins to 
fall, w represents the angular velocity of the earth's rotation, <^ the latitude 
of the point Po , W the potential function of the statical field of force which 
is at rest with respect to the rotating earth, and the symbols W\ , W^^ , W\^ , 
W\( , W^^ , W\^ , stand for the values of the derivatives 

dW d^W d^W d^W d^W d^W 

at the point Po . If we put — ^ = h = RPo , where h is the height through 
which the body falls, the corresponding values of 17 and | are respectively the 
easterly and the southerly deviations which correspond to h . 

Let us now express ^ and $ explicitly in terms of h. The last equation 
assumes the form 

a 
where 

a ^W\, y = ^{W\i - 4.0^ cos" <i>) , 6 = - ^ co cos<^ • IF^^ . 

Hence 



i 



a 



and this when solved for t yields the relation 



+ 



Since t is positive, the upper sign must be used. The expressions for 17 and | 
then assume the forms: 

, = 3T/2 • cocos.^ . ^^==+g ^— j^ A^ 

l/2 
+ ^ CO (2 sin <^ • r|f + 2 cos <^ • TF;, - 3 cos <^ ■ W\i 

15/2 

II. +4«-(5cos^<^-2cos2<^)) 



(- W\f>^' 

1 F 

^ = - ( 4a)^ sin (> cos (> + 51F^^) / _ p^o x 



+ -jy w ( 9 sin <^ IF^f + cos <^ • F?, ) T^^o yn ' 
and the expression for f assumes the form f = — ^, as it should. Since 
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JF, , W^, W^ are the rectangular components of the acceleration g due to 
weight, we may write 

-.■•=,. ,n<=-(^').' """'-(Ml '^«-(0).- 

and therefore the above formulas may be written as follows : 
, = 3l/2.cocos<A-^,-g 1 



V2 
30 



+ -43cos*(0\-2,in*(ffX+2eo.*(^)^ 



III. + 4co2 ( 5 cos' <^ - 2 cos^ ) 

c 4to^ sin cos .^ - 5{dg/d^)o 






+ ll4^^^*(^)o-^^^'^K'^')o] 



^2 



For a distribution of revolution we have 



= 0, T-r.= 0, 



and therefore the term in I? in the expression (rt) for the easterly deviation 
and the term in A^'^ in expression (| ) for the southerly deviation drop out. 

The establishment of the general equations I. 

§ 2. The curve c. The curve c has already been defined as the path of the 
falling particle with respect to the rotating earth. Let us first refer it to a set of 
rectangular axes 0—x,y,z, fixed in the rotating earth, and such that the origin 
is a fixed point (interior to the earth) of the earth's axis a of rotation. Oz 
is coincident with a, and positive in the direction of the north pole, Ox is 
perpendicular to Oz and is so chosen that the initial point Po lies in the plane 
zOx, the positive direction being that from Oz to Po, and finally Oy is per- 
pendicular to Oz and Ox and is positive in the direction which Ox would have 
if it were revolved around the axis Oz through an angle of 90° in the direc- 
tion of the earth's rotation. (See Fig. 1.) 

If we represent by co the angular velocity of the earth's rotation, and by V 
the function f{x, y, z) which is defined by the integral : 

U= f~, p=Vix- XiY +{y- Vif +{z - Zif, 

oy P 

where dm represents that element of mass of the earth which is situated at 
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the point whose coordinates, with respect to the axes — x,y ,z, are a;,- , yi , Zi , 
the integration being extended throughout the whole volume of the earth, 
then the function: 



(1) 



W= f7+2-(^'+y') 



is the potential function of the statical field of force which is at rest with respect 
to the rotating earth* (i. e., the lines of force: 



dz : dy :dz = 



dW dW dW 



dx ' dy ' dz 
of this field are at rest with respect to the rotating earth f)- 




Fig. 1. 

In order to find the differential equations of motion, let us denote by 
— a;i , t/i , 2 a set of rectangular axes which are fixed in space and so oriented 
that at the instant t = 0, when a falling particle leaves Po (at rest with respect 
to the earth), the set of axes — x, y , z coincides with the set 0—Xi, yi, z 
(see Fig. 1). If we denote by d the angle XiOx, positive in the direction in 
which the earth rotates, the relations between the coordinates are: 



Then 



X = xi cos -\- yisin d , y = — xi sin 6 -^ yi cos d , z = Zi . 
U = f{x,y,z) = fi{xi,yi,z,t), 



* See PizzETTi, TraMato di Geodesia teoretica, § 2 (1905). 

t It should be borne in mind that the level (equipotential) surfaces, W = const., are not 
necessarily surfaces of revolution. 
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where 6 = wHs a function of the time t . The differential equations of motion 
with respect to the fixed axes — Xi, yi, z, are 



d^xi dU 


d'y^ dU 


d^z dU 


d^ dxi ' 


df dyi ' 


df dz ■ 



If we subject these differential equations to the inverse of the transformation 
just given, namely to the transformation 

xi = X cos 6 — ysin d , yi = x sin -\- y cos d, zi = 2 , 

we obtain a set of differential equations* which is equivalent to the set: 

d^x ^ dy , dU d^y , ^ dx , dU dH dU 

Since by equation (1), 

aZ7 , dW dU , dW dU dW 

these equations assume the forms: 

^■' df '^dt~~d^' df'^dt~dy' ~df~dz- 

The curve c is then that solution of the differential equations (2) which is 

subject to the initial conditions: 

/ON TTTi <^ ^ dx dy dz ^ 

(3) When t = 0, x=ro, y=0, z = Zo, 'dt^H^dt^ 

§ 3. The curve c referred to the cardinal axes at Po ■ Before we can define 
the cardinal axes it will be necessary to define the following terms. At a 
general point P which is at rest with respect to the rotating earth, the vertical 
is the straight line which coincides with the plumb-line,t the plumb-bob of 

* <Px 

dU a dV . ^ 

— -i~ cos B r— Sin 6, 

dx dy ' 



- cos 9 — T^ sin e — 2 ( ^ sin e + ^ cos fl j 0) — ( a; cos e — y sin e ) co2 



(Px . 



e + -Tf COS e + 2{ -T:cose — -^sine] a— (a;sme + «cos9)a)'=— -sin9 + — — cosff. 
aP \ at at J dx ay 

If we multiply the first of these equations by cos and the second by sin and add, we 
obtain 

d^x dy dU 

and if we multiplj' the first by — sin B and the second by coa 6 and add, we obtain 

d^y,^ dx . , dU 

t It is assumed that the plumb-bob is a heavy particle and that the line is weightless and 
perfectly flexible. 
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which is situated at P (this line does not necessarily intersect the axis of 
rotation a of the earth), the meridian -plane \s the plane which passes through 
the vertical and is parallel to the axis a , the horizontal plane is the plane which 
passes through P and is perpendicular to the vertical, the east-and-west line 
is the straight line which passes through P and is perpendicular to the meridian 
plane, the north-and-south line is the line of intersection of the meridian and 
the horizontal planes. By the cardinal axes at the point P we shall mean the 
vertical, the east-and-west and the north-and-south lines at the point P. 




Fig. 2. 



The positive directions of these axes shall be those toward the zenith, the east, 
and the south respectively. A few more terms will now be defined. The 
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latitude at P is the complement of the angle which the positive direction of the 
vertical at P makes with the positive direction of the axis a . The longitude 
at P is the angle which the meridian plane at P makes with a fixed plane through 
the axis a (we shall here assume that the fixed plane is the plane zOx ) measured 
from 0° to 360° to the east.* 

Let us now take the cardinal axes at the point Po as the axes of the co- 
ordinates ^, r), f, where ^ is measured along the north-and-south line, t] 
along the east-and-west line and f along the vertical. See Fig. 2, in which 
Oz = a represents the axes of rotation of the earth, and Ox, Oy & pair of 
t. 




Fig. 2'. 



rectangular axes which are at rest with respect to the rotating earth and are 
perpendicular to Oz at . Pq is a point of the plane xOz . Po M is parallel 
to a , and intersects the plane xOy in M . We put 



'•o 



OM, Zo=MPo. 



The Hne Po L is the vertical at Po ; it intersects the plane xOy in L . LPo M is then 
the meridian plane at Pq . The plane NESW , which passes through Po and is 
perpendicular to Poi , is the horizontal plane at Po . These two planes intersect 
in the north-and-south-line NS . Let us denote by t the plane which passes 
through Po and is aparallel to the plane xOy . This plane intersects the hori- 
zontal plane in the east-and-west line EW . 

We agreed to take the plane xOz as the plane from which t'^ measure longi- 
tude. Therefore the longitude of Po is the angle AX = ;^ OML . The latitude 
of Po is the angle <^ = i^ MLPo . 

* See PizzBTTi, loc. eit., § 5. 
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In order to get the equations of transformation from the axes — x, y , z 
to the axes Pq— f , t; , f , let us introduce the auxiliary variables €i , f i , ^2 , 172 , 
where |i represents distance measured along ML from M in the direction LM , 
and f 1 , f 2 and ??2 are defined by the equations 

Then 

?2 = li COS AX — 17 sin AX , 172 = ^1 sin AX + j; cos AX , 
and 

^1 = f cos <^ + ? sin <^ , f 1 = f sin <^ — ? cos (^ . 
Therefore 

X = ro + ( f cos (j) -jr ^ sin <^ ) cos AX — 17 sin AX , 

(4) y = ( f cos <^ + ? sin <^ ) sin AX + rj cos AX , 

2 = Zo + f sin <^ — ^ cos <l> . 

The inverse of this transformation is the transformation: 

I = [{x — To) cos AX + y sin AX ] sin <^ — ( 2 — Zq ) cos <^ , 

7} = — (x — ro ) sin AX + y cos AX , 

f = [ ( X — ro ) cos AX + 2/ sin AX ] cos <^ + ( 2; — 2o ) sin <^ . 

If now we subject the system of differential equations (2) to the transformation 

(4), we obtain the system of differential equations,* 

*The details of this process are the following: 

^_ 2 ^ =^^ 
dp dt dx 

becomes 

cos AX I ~ cos <^ + -3i| sin (^ j — sin AX -3^ — 2 sin AX ( -3^ cos ^ + ^ sin (^ j w — 2 cos AX • t^ . « 

dW . ^ ., dW . ,^ , dW 

= -TT- sin d> cos AX r — sin AX + -tt" cos 4> cos AX. 

3J drj df 

dp "*■ dt dy 
becomes 



Ax(Jfcos<^ + ^sin^) + cosAx|^+2cosAx(f cos^+gsin*)«-2sinAxg, 



■ -iTT^ sin <^ sin AX + — — cos AX + -^r '^°^ "^ sm AX . 
o« . "'7 OS 



sin 



becomes 

sm <j> . 

If we multiply these equations respectively by cos AX sin 41, sin AX sin , — cos ^ and add^ 
we obtain the first of the equations (6) ; if we multiply respectively by cos A X cos <j> , sin AX 
cos <i> ,sm<l> , and add, we obtain the second of the equations (5) ; and if we multiply respectively 
by — sin AX , cos AX , , we obtain the last of the equations (5). 





O'z dW 
dP dz 




^^,COS*= ^^COS0+^^ 
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(5) ^+2a,(^sin<>-^+cos.^-^j = 



dW 
dt] 



(P^ „ dv dW 

-^-2a,cos<^.^=-^. 

The initial conditions (3) assume the form: 

(6)me„ ,-0, «=,.f-0. S4' = |-0. 

Therefore the curve c is that solution of the differential equations (5) which is 
subject to the initial conditions (6). 

§4. Another derivation of equations (5). The differential equations (5) 
may also be obtained as follows. The projections on the moving axes i? , ? , f 
of the absolute velocity Va are * 

Vai=J^-VVl+rn-v^, 

where F°, V\, V\ are the projections on the axes tj, ^, f of the absolute 
velocity F" of the initial point Po , and p, q,r are the projections on the same 
axes of the rotation &> of the earth. From Fig. 2' it is easily seen that 

r° = F" COS AX, ¥1= F« sin AX sin 4,, V} = ¥*> sin AX cos .^, 
where 

F« = ccro, 
and 

p = , q = 01 cos 4>, r = — (a sin <t> . 
Therefore 

Var, = jX + "''0 COS AX + CO COS (^ • f + CO Sin (^ • ^ , 

(7) Foj ~ liJ ~^ '^''"'^ ^^^ AX sin <^ — CO sin (^ • 7j , 

d^ , 
Vai = "jT + '^'"o sm AX cos <^ — co cos <t> • v- 



'See Appell, Traiti de MScanigue Raiionnelle, vol. 1, (1902) §61. 
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The projections on the moving axes ?; , ? , f of the absolute acceleration J a are 

Jar, = -^ + qVai " tV a^ , 
Ja^ = ^-^ + rVa„ - pVai, 
Jai = -^+ pVai- qVar,. 

Hence by the expressions for Fo, , Va( , Va^ , 

Ja,, = -^+ 2cof cos<^ •-^+ siiKt,-^]- u^ (v - rosin AX), 

(8) Jaj = -^ — 2w sin <^ ■ -T7 — 0)2 sin (^ (cos (^ • f + sin <^ • 1+ ro cos AX), 

Jai = -jw — 2ia COS <t> • -^ — (0^ cos <t>{cos<t> • f + sin <> • ? + ro cos AX) . 

Since gravitation acts instantaneously, the projections of the absolute force 

on the axes »?, ^, f are 

dU dJJ dJl 

dv' 3| ' af • 



(9) Jo, — -,_ , Ja^ ~ at ' "f ~ ai- ' 



Hence the differential equations of motion are 

dU _dU _d_U 

dv ' "^ ~ d^ ' ^"^ af 

where the expressions for Jo, , Joj , Ja^ are given by (8). In virtue of relations 
(1) and (4) we find that 

— — = — - + of iv— rosmAX) , 
arj or] 

-J Tj/ a TJ 

— — = —- + co^ sin <^ (cos <^ • f + sin <^ • ? + ro cos AX) , 

o TXT' fj TJ 

-J— = — - + 0)2 cos <^( COS <^ • f + sin <^ • ^ + ro cos AX) , 

and hence equations (9) are identical with equations (5). 

§ 5. Integration of the differential equations of motion (5). Let us write 
equations (5) in the following form : 

r," + 2co (sin <t> ■ ^' + cos4> ■ t') - W^= , 

(a) ^" - 2o}sm(t> ■ v' - W( = 0, 

f " - 2o) cos <t> • v' - Wi = 0, 
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where the primes (') and seconds (") denote the first and second derivatives, 
with respect to the time t, of the functions to which they are attached, and 
the subscripts (,,{,;) denote the first partial derivatives of W with respect 
to the attached subscript. In the equations which follow thirds ('"), etc., and 
double subscript (,,,,{ . . .), etc., will be used to denote the higher derivatives. 
For the particular values of i; , ^ , f , and their derivatives with respect to t , 
which correspond to < = , we shall use the symbols tjo , ?o , f o > Vo> ^'o' f o > 
Vo ) C ' 6tc. Hence condition (6) may be written in the following form : 

For the particular values of the derivatives of W which correspond to t? = jjo , 
I = ^0, f = fo, we shall use the symbols PF», W^, W\, W^^, W^^, W^^, 
Wh, W\^, W\i, etc. 
For a set of functions 

r,= \{t), ^=m(0, ^=v{t), 

which is a solution of the set of differential equations (a ), the left hand member 
of each of the equations (a) vanishes identically. Therefore, for this solution, 
the derivative with respect to t of each of these left-hand members must vanish 
identically. Hence the identities: 

„"' + 2co(sin,^ ■ r' + cos,^ ■ r)-[W^ ■ 7,'+r,f ■ r+TF,f • f']=0, 

(a') r"-2cosin<^-V' -[W^f,-V+W^«-r+W^«-f']=0, 

r'-2cocos.^- V -[W"f, -V+W^a- r+TFff •f']=0. 

Similarly, if we assume the existence and continuity of all the derivatives 
which are needed, we obtain the following sets of identities : 

n'" + 2co (sin 4, • r" + cos <^ • f") 

- [ ^„ • v + w^^ ■ r' + w,^- f" + ai V + ^1 r + 71 n = o, 

(a") ?'^ - 2co sin <^ • tj'" 

- [ TFj, • n" + w^« • ^" + w^^ ■ f" + a2 V + ^2 r + T2 n = 0, 

f '^ - 2co cos <i> • ri'" 

- [W^, ■ v" + W^i ■ ^" + W^i ■ r + a,v' + Ps r + 73 f] = 0, 

where «», j3,- , yt, i= 1, 2, 3, stand for homogeneous linear expressions in 
v' > ?' > f ' ) the coefficients of which are third derivatives of W* , 

*J(u—u{x, y , z) , u' = Ux • x' + Uy • y' + Ui z' , 
u" = Ux ■ x" + Uy ■ y" -\- Uz • z" -\- { Uxx x' + Uxy y' + Ux. z' ) x' 

+ ( Uyx X' + Uyy J/' + Uy z «' ) ?/' + ( U zx X' + U zy V' + U,, Z' ) Z' 
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rf + 2w (sin <^ • ^'^ + cos <(> • f '^) 

- [W^ ■ v'" + W,^ ■ r + W,i ■ r'" + a[v' + P[ r + 7'i f] = 0, 
(a'") ^ — 2wsin<^ • 17'^ 

-iw^,- v"'+ jr« • r'+ Wii ■ r+a, -v'+p', r+ 7; f] = 0, 

r — 2wcos<^ • j?'^ 

- [w^,- v"'+ w^« • r+ w^i ■ f"'+«; v'+P's r+ 7; r'] = 0, 

where «• , /3j , 7', , i = 1,2,3, stand for polynomials in the »?' , I' , f' , v" > ?" > f " 
and the third and fourth derivatives of W.* 
By condition (/3) we already know that 

'?o = ?o = fo = , i/o ~ ^0 ~ f ~ • 

We also know that 

U7»= JF| = 0. 

Therefore, when t = 0, the identities (a) yield the relations 

'?o=0, ^„' = 0, ro = TF|; 

the identities (a'), the relations 

V;= -2cocos<^-IF^ c = o, r;" = o; 

the identities (a"), the relations 

i^ = ^? Ki > io = fT? ( Wli - ice' sin <#, cos .^ ) , 

f''= PFnW^?f-4co2cos2<^); 

and the identities (a'"), the relations 

,; = - 2co Wl [sin4> ■ Wli+ cos <t>iWli+Wl - 4o,')], 

?J = + 2co TF?[sin<^ • Wli - C0S4,- Wl,], 

n= + 2co PF^cos<^[ TF;^ - Wl] = 0. 
and 
u'" = M;r x'" + uy y"' + u, z'" + 3 [ ( M„ x" + u,y 2/" + M„ 2" ) x' 

+ ( Mxl; X" + Myj, 2/" + M„. 2" ) 2/' + ( Mji l" + Mj, 2/" + Mj, 2" ) «' ] 

(MxxxO:' +M„„ 2/' +«»»««' )«' "] r («»»»«' +%w2/' +"»i»2')a:'1 

+ -{ +(Ma^Bia;'+Mii,i,2/' + M*»«2')2/' f^'^' + T +(Mraia;' + M„„„^'+M„„.2')2/' Vl/ 

+ (m«2xX' +Mi«v2/' +Wx«z2')2' J L + (%.xX' + M„,„2/' +M„„2')2' J 

(w«lX' +Uxyty' +U!x.z')x' ^ 

+ ■{ +iui„xx' +u,y,y' + u.tiz')y' >z'. 
+ (uiixx' +U::yy' + Uz!.z')z' J 
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Let us now assume that the conditions are satisfied under which the set of 
differential equations (5) or (a) have a solution of the form 

ri=ao+ait+a2fi+a3f+ait^+aifi+ ■■ ■ , 

(10) ^=bo+ht+b2f+hf+bif+b,f+---, 

^ = Co + Cit + Cif + Csf + Cii* + c,f + ■ ■ ■ , 

in the neighborhood of the point Po. It then follows, from the preceding 
work, that for the initial conditions (6) or (/3), 

ao = 170 = 0, «! = 179 = 0, a2 =J v'o' = , 
6o=lo=0, 6i=^; = 0, 62 = i|;'=0, 

co = ro = o, ci=f; = o, c2 = Uo' = Wh 

as = Wo" = — if} COS <t> • Wi , 
bz = Wo" = 0, 
C3 = Wo" = ; 

hi = irV^i' = ^hW} ( Wli - 4co2 sin 4> cos ^) , 

C4 = 2Vf«' = ifW'i ( W'^i - 4«2 cos'' ./.); 

«5 = lio 1?; = - eVwW^f [ sin « • W^?f + cos <^ ( JT^j + W',, -ic')], 

h = ik^o = eV'ol^? [ sin <^ • PF«j - cos <^ JF|,] , 

cs = T^-jrfo — • 

The equations (10), for which the constants are given by the relations (11), 
are then the equations of the curve c referred to the cardinal axes at Po • 

§ 6. The plumb-bob locus referred to the cardinal axes of Pq . This curve has 
already been defined as the locus of plumb-bobs of all plumb-lines which are 
supported at the initial point Po, of the falling particle. It is therefore the 
locus of the feet of perpendiculars dropped from Po(i7 = 0, 1 = 0, f = 0) 
to the level (equipotential) surfaces W = const. Its equations are easily seen 
to be 

V _ ^ _ ^ 



(11) 



<12) 



dW dW dW 



dr, a? df 

Now let us solve these equations in the form 

17 = ai f -1- a2 r' + • • • , ? = ;8l f + 182 f' + 

Trans. Am. Math. Soc. 31 
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In order to do this let us write the equations (12) in the form: 

Fi(v, i,r) = o, 



(i = l,2) 



where 



^1 = 






and 



F.= 



W, 



V 



The derivatives dri/d^, d^/d^, d^ri/ d^^, d^ ij d^'^ are expressible in terms 
of the derivatives of i^,- , i = 1 , 2 ,* and the derivatives of Fi are expressible 
in terms of the derivatives of W .f The coeflBcients ai, a^, /3i , /32 have then 
the following values: 



ot\ 






wy 



02 = 



wy 



Therefore the plumb-bob locus may be represented by the equations 
(13) 



r. = — ^ f2 I 
'? TJ/O S , ^ 



r? 



? Tl/o S T 



TF? 



§ 7. The curve c referred to axes of origin Po and directions which are parallel 
to the cardinal directions at the plumb-bob of a plumb-line which is supported 
at Po . Let us denote by R (Fig. 3) the point at which the plumb-bob of a 



* The relations are the following: 



dv df ^ 3| df ^ 3f 



dFi (Py 
dv df'' 



.dr^r d'Fj dl d^ dji 
d^ldvdi df"*" dv'' df 



+ 



d^Fi 






difd^di^ d^dy d^i ,a^di ^ 
dr L de ar d|a, df ■'"afad '^d^did^ 



t The relations which are needed are the following: 



dFi 
dv 




r 

Win' 


dFi 


1 f 

Wu Wis 


+ 


1 

Wi Wi 


aFi 


1 r 


dFi 
dv 


r 

Win 


Wr„ + 


1 dFi ^ V dFi 

Wi Wr, ' di Wi$ Wr,( ' 3r 


r 1 

Wii Wni 






d^F, 1 r 

dt- WHi Wiii 


,01 „ 

"*" Wh Wii ~"^^^' 








d^ 


Wiii w 


iCC 


+ 1 


T^.. +^'' 


t 



-t- 



(i = l,2), 



(i = l, 2). 



I ° ^ 

Itfj Wi 

I 1 



Tf f Wr, 



(f).=«. (f).= <. (f).= -'^-». 



3f« 



= -2Tr|i, 



%■'-<<■ 
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plumb-line which is supported at Po is in equilibrium under the tension of 
the plumb-line and the force of the field whose potential function is represented 
by W . The line RPq is then the vertical at R . Let us denote by RE jt and RSg 
the easterly and the southerly directions at R . Now let us draw through Po 
the axes Pov, Pol, Pol parallel to REg, RSg and RPq respectively ( Pq ? 
is the continuation of the plumb-line PPo) . Let us denote by (ai, /3i, 71) , 
(«2, P2, 72) and (az, p3, 73) the directional angles of the axes Po — V > 1,1 
with respect to the cardinal axes Po — 77, ^ , f . Then 



(14) 



17 = cos ai • T) -\- COS /3i • ^ + cos 7i • f , 

1 = COSai • 1)+ COS /32 • I + COS 72 • f , 
f = COS as • 77 + COS /33 • I + COS 73 • f . 




Fig. 3. 
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In order to obtain expressions for these directional cosines let us recall that, 
at a general point P , (1) the vertical is the normal to the level (equipotential) 
surface which passes through P , (2) the east-and-west line is the intersection 
of the horizontal plane and the plane (z = const.) which passes through P 
and is perpendicular to the axis of rotation of the earth, (3) the north-and-south 
line is the common perpendicular to the vertical and the east-and-west line. 
Therefore, at a general point P, (ri, f, f), (1) the directional cosines* of 
the vertical are proportional to 

w,, w^, WiV 

(2) those of the east-and-west line are proportional to 



where 



H = 



w, 


w, 


dz 


dz 


3? 


dt 



H, S, Z, 



sin <j> • W(-{- cos<j) ■ W^, 



w, 


w. 


dz 


dz 


ar 


dr, 



— sijKj} ■ W. , 



Z = 



w. 


W, 


dz 
dti 


dz 
d'k 



— cos <t> • Wr,; 



and (3) those of the north-and-south line are proportional to 



w, 


W, 




w, 


w. 




w. 


W, 


E 


z 


* 


z 


H 


9 


H 


s 



If the point P is not too far from Po , we may write without appreciable error. 



Therefore, for the point R , whose coordinates are 






wi, 



w 



_Ek,. 



'' n/9 ^ ' ^ IT/0 r > f r ; 



TF? 



we may write 



W^ = wii 



f, 



wf = w; 



a 



f, Wf=W\+W\i-i, 



' With respect to the cardmal directions of Po . 
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and hence 

cosa3 :cos^3 :cos73 = PF»f • f : W^^ • f : W^H ^^ff " f' 
cos ai : cos /3i : cos 71 

= cos(t>-W'i+ {sm<l>-Wli+ cos <I>-Wli) r.- sin<t,-W'^i- r-- cos<t> -W'^i- ^ , 

cos a2 : cos 02 '• cos 72 

= sin.^ • PF? • PF,»f • f : cos.^ (JF^)^ 

+ JF^ (sin • r^f + 2 cos .^ • Wii) f : - cos <^ • JF? • TF|f • f, 

for terms of order not greater than the first. In order to get the cosines them- 
selves, we observe that * 

((n • f)^+(n-f)^+(^?+nr)^)-* = 1^-0(1 -^'f +•••), 

({cos..^- TF»f+(sin,^. W]i+cos4> ■ TF|f ) f p + (- sin<^ • W^^ ■ f)^ 

+ (-cos<^.pF;,.r)^)-*=^^^— f^„(i cosVtf^ — ^^^+--7' 

( (sin .^ • W'( ■W'r,i- n'+{oos<j>- { W'^r + Wl (sin <l> ■ Wj^ 

+ 2cosct>- TF^f)r}2 + (-cos.^- w'i-wii- n')-* 

sm<t> ■ W'(i+2 cos </. • W^i 

cos4> ■ {W}yV 

Therefore 



I ( 1 _ '<' " " ' -^ ll-4> J- -I- 

cos<t>-(W}y\ cosct>-Wl ^^ 



cos as = "s^ f , cos 03 = -^ f , cos 73 = 1 , 

We W\ 

(15) cosai=l, cos /3i = — tan (^ -^ f , cos 71 = — ^j^ • f , 

cos a2 = tan <^ • -jj^ • f , cos /32 = 1 , cos 72 = — -frMT • f > 

for terms of order not higher than the first. In virtue of the values (15), 
equations (14) assume the form: 

^ = 5j — tan (^ -^ ?f - -^ fS 

TF% PF2^ 

(16) 1 = tan.^.-j^lSf+?--^rS 



f- pp« ''^"'"ir? ^^''' ^' 



f F-K f 



T7T Ai\ 2A' 



f+— ?7^ f'+' 



i/A + Bf + Cf2 + ■ • • ^H 2 A ' 8^2 
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in which we must now substitute for )? , ? , f the values in t which are given by 
equations (10). Doing this, we find * 

n=-^o>cos<t>-W'i-f-^j Wl • TFJf • ^ - ^\ CO W^i (sin .^ • TF|f 

(17) -i= - (i cousin ,!> cos <!> + ^^■^Wli)W} ■ i* 

- ^\ CO W'i{9 sin <j>Wli+ cos <t> WD f, 

f = ^Wl ■ i' + j\ (Wli - 4 co^ cos" <l>) W} ■ i^ - ^ CO cos <t>W}W'r,(f. 

Equations (17) are the equations of the curve c referred to the axes 

If we put — f = h = RPo, where h is the height through which the particle 

falls, then the corresponding values of rj and | are the easterly and the southerly 

deviations which correspond to h. These, and not the quantities t;, ^, are 

the ones which are measured in experiments. 

Washington Univbesitt, 

St. Louis, Mo., A-pHl,. 1912. 

* We can not carry these developments beyond i' unless we carry the developments (15) 
beyond f. 



